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Abstract 

In this letter we present the analytic results for the two-loop corrections to the Wilson 
coefficients Cg(fiw) an d Cio(/ i vy) i n type- 1 and type-II two-Higgs-doublet models at 
the matching scale nw These corrections are important ingredients for next-to-next- 
to-leading logarithmic predictions of various observables related to the decays B — > 
X s l + l~ in these models. In scenarios with moderate values of tan/3 neutral Higgs 
boson contributions can be safely neglected for I = e, \i. Therefore we concentrate on 
the contributions mediated by charged Higgs bosons. 



1. Introduction 



In the standard model (SM) rare decays of B-mesons like B — > X St( f) or B — > 
X s ,dl + l~ are induced by one-loop diagrams. In many extensions of the SM, there are 
additional one-loop contributions in which non-SM particles propagate in the loop. If 
the new particles are not considerably heavier than those of the SM, the new contri- 
butions to these decays can be as large as the SM ones. As an illustration of the high 
sensitivity of these decays to new physics, we mention that the most stringent bound 
on the mass of the charged Higgs-boson in the type-II two-Higgs-doublet model comes 
from rare 5-decays, viz B -> X s7 , leading to M H > 280 GeV (99% C.L.) [1]. 

It goes without saying that one should try to get information on the parameters in a 
given extension - the two-Higgs-doublet models in this letter - from all processes which 
allow both a clean theoretical prediction and an accurate measurement. This means 
that precision studies similar to those for B — > X s ^ [2-5], where higher order QCD 
corrections are crucial, should also be done for the process B — > X s l + l~ . On the the- 
oretical side this means that next-to-next-to leading logarithmic (NNLL) calculations 
for the branching ratio and/or the forward-backward asymmetry are needed. 

In this letter we consider QCD corrections to the process B — > X s l + l~ (/ = e, fi) 
in 2HDMs. We neglect diagrams with neutral Higgs-boson exchange. This omis- 
sion is justified in the type-II model, if the coupling parameters (mi/Mw) tan (5 and 
mi/ \M W cos (3) are sufficiently smaller than one. In this case the operator basis is 
the same as in the SM. Only the matching calculation for the Wilson coefficients gets 
changed by adding the contributions where the flavor transition is mediated by the 
exchange of the physical charged Higgs boson. While these extra pieces are known for 
the coefficients C-j and C% to two-loop precision for quite some time, the correspond- 
ing results for C 9 and C w , presented in this letter, were not published before. The 
phenomenological consequences for the branching ratio and other observables will be 
discussed in [6]. 

The remainder of this letter is organized as follows: In section 2, we summarize the 
necessary aspects of the 2HDMs. In section 3 we first present the effective Hamiltonian, 
followed by the analytic results for the charged Higgs boson contributions to the Wilson 
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coefficients Cg(/iw) and Cio(a*w)- In this section we also briefly investigate how the two- 
loop corrections reduce the renormalization scheme dependence related to the definition 
of the top-quark mass. 

2. Two-Higgs-doublet models 

In the following we consider models with two complex Higgs-doublets 0i an d 02- 
After spontaneous symmetry breaking these two doublets give rise to two charged (H^) 
and three neutral (H°, h°, A ) Higgs-bosons. When requiring the absence of flavour 
changing neutral currents at the tree-level, as we do in this paper, one obtains two 
possibilites, the type-I and the type-II 2HDM [7]. The part of the Lagrangian relevant 
for our calculation is the Yukawa interaction between the charged physical Higgs bosons 
and the quarks (in its mass eigenstate basis): 

c ^M(jt) XWL ' v '- d «- + Ga) y w^K +h - c - (1) 

The couplings X and Y are 

X = -cot/3, Y = cot/3 (type-I), 
X= tan/3, Y = cot/3 (type-II), 

where tan/3 = i^/fi, with v\ and v 2 being the vacuum expectation values of the Higgs 
doublets 0i and 02, respectively 

In the following we will use the generic form (1) for the interaction between H ± 
and the quarks. It will turn out that the Wilson coefficients Cg(nw) arid Cio(f^w) are 
independent of the model (type-I or type-II), as they only depend on Y 2 . 

3. Charged Higgs contributions to Cq(^lw) and Cio(^w) 
at the two-loop level 

In this section we first briefly describe the effective Hamiltonian. We then present 
the analytic results up to two loops for the charged Higgs boson contributions to 
Cg(fJ>w) and C w (fiw)- Finally we briefly investigate the impact of the new two- loop 
contributions on Cg([/,w). 
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3.1. Effective Hamiltonian 



To describe decays like B — > XJ + l~ we use the framework of an effective low-energy 
theory with five quarks, obtained by integrating out the heavy degrees of freedom. In 
the present case these are the t-quark, the and Z° boson as well as the charged 
Higgs bosons H ± , whose masses M H are assumed to be of the same order of magnitude 
as M w . As in the SM calculations we only take into account operators up to dimension 
six and set m s = 0. In these approximations the effective Hamiltonian relevant for our 
application (with \AB\ = \AS\ = 1) 

AC 10 

n eff = ~^0*E c iWW ( 2 ) 

contains precisely the same operators Oi(fj) as in the SM case. They read: 
0i = {s Ll ,T a c L ) (c L ^T a b L ) , 2 = (s Llll c L ) (c l1 %l) , 

3 = (s Lllx b L )E q (qrq), 4 = (s Llfl T a b L )E q (qrT a q), 

5 = (s Ll ^ ulp b L )E q (qrYYq), O e = (s Ll , lvlp T a b L )Y, q (qYl u l p T a q), 
7 = j>m b (s L <j^b R )F^ : O s = j-m b (s L a^T a b R )G^, 

O g = $(s L ^b L )Ei(iri), O 10 = 4(^7^)^(^750, 

ij s y s 

(3) 

where T a (a = 1, ...,8) are the SU(3) colour generators, and g s and e are the strong 
and electromagnetic coupling constants, q and / appearing in the sums run over the 
light quarks (q = u, ...,&) and the charged leptons, respectively. 

The Wilson coefficients Cj(/i) are found in the matching procedure by requiring 
that conveniently chosen Green's functions or on-shell matrix elements are equal when 
calculated in the effective theory and in the underlying full theory up to O [(external 
momenta and light masses) 2 /M 2 ], where M denotes one of the heavy masses like Mw 
or M H . The matching scale /iw is usually chosen to be at the order of M, because at 
this scale the matrix elements or Green's functions of the effective operators pick up the 
same large logarithms as the corresponding quantities in the full theory. Consequently, 
the Wilson coefficients Ci(fi w ) only pick up "small" QCD corrections, which can be 
calculated in fixed order perturbation theory. For the following it is convenient to 
expand the Wilson Ci(nw) as 

Ci(nw) = Cf W) + (J^W) + (J^ (2 W) + 0{gl). (4) 



We note that due to the particular convention concerning the powers of the strong 
coupling constant g s in the definition of our operators, the contributions of order g 2n 
to each Wilson coefficient originate from n-loop diagrams. 

In the SM all the Wilson coefficients Ci(/j, w ) are known at the two-loop level. In 
2HDMs, the charged Higgs boson exchanges lead to additional contributions. For the 
following discussion, we split the Wilson coefficients into a SM- and charged Higgs 
boson contribution according to 

C.(fi w ) = C i:SM (vw) + Cifl(ti w ) . (5) 

The individual pieces Ci^mi^w) and Ci^il^w) can be expanded in g s in the same way 
as Ci(fj,w) in ec i- (4). While C^^il^w) and Cs,n{f^w) are known at the two-loop level 
[2,3,5], Cg^u(^w) and Cw^l{^w) were up to now only known to one-loop precision [8,9]. 

(2) (2) 

3.2. Analytic results for Cg^fiw) and C) ^ h (/xvk) 

(2) (2) 

We did the matching calculation for Cg ^fiw) and C{ ^ u (fiw) in two different ways, 
leading to identical final results: On the one hand we performed a matching calculation 
for (the off-shell) Green's function related to b — > sl + l~ , as described in detail for the 
SM in [10]. On the other hand we matched the corresponding on-shell amplitude onto 
the effective theory, following basically the methods described in [11], but using some 
simplifications*. In both methods, the hard part of the calculation consists of working 
out the one-particle irreducible diagrams shown in fig. 1. 

After using heavy mass expansion techniques [12], partial fraction decomposition 
and the usual reduction of tensor integrals to scalar ones, we obtain integrals of the 
type 

(2) ( m 2 ) n 1 +n 2 +n 3 -4+2e . ^4-2^ ^4-2^ 

nin2 ™ 3 " tt 4 ^ T(l + e) 2 J (q\ - m 2 )"i {q\ - m 2 2 ) n * - q 2 ) 2 ]^ ' 1 ' 
which are known explicitly [13,14]. 

Note that only the contributions from the internal top-quarks have to be taken 
into account in these diagrams, because the charm contributions, which come with a 



*As a byproduct of our calculation, we also confirmed the known result for the charged Higgs 
contribution C) ^ (see e.g. [2,3]). 
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Figure 1: One-particle irreducible two-loop diagrams for b — > sl + l needed to extract 
the charged Higgs boson contribution to Cg (f^w) and C[ (f^w)- The external quark 
lines (solid) denote the incoming 6-quark and the outgoing s-quark, while the wavy line 
denotes a virtual photon or a Z°-boson, which decays into a / + Z~-pair. The internal 
dashed-, solid- and curly lines denote the charged Higgs boson H ± , the t-quark and 
the gluon, respectively. 

relative suppression factor of mbm c /mf or ml/mf, only induce dimension 8 operators 
which are neglected in our treatment. 

We write the one- and two-loop charged Higgs induced contributions to Cg(/iw) an d 
C w (/iw) m the form 

d&M = y2 ( r(n) + H > 

V s w J 

= -Y^(^Z^), (7) 
V s w ) 

where s w = sm9 w . The terms proportional to (T^) account for the n-loop Z°- 
(photon-) penguin diagrams. 
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The one-loop contributions and read [9] 

r( i) -(38 - 79y + A7y 2 )y (4 - 6 y + 3 y 3 ) y } 

108 (y-lf 18(y-l) 4 

Z (D = xy ^y\ny, (8) 

8(y-l) 8(y-l) 2 i; 



with 



2 2 



X ~~ ' f 71^2 " (9) 



Note that and depend via x and y on the renormalization scheme for the 
t-quark mass. To illustrate this dependence, we give our results in the commonly used 
MS- and pole-mass scheme. The relation between these mass definitions is given by 

— / \ pole A . 2a 5 (^ w ) m]? olc Aa s (fi w )\ 2 
m t (ii w )=m i t 1 + In )+0(a s ), (10) 

\ 7T ill) 

where rn t {^w) and m^ olc are the top-quark mass in the MS-scheme and pole-mass 
scheme, respectively. 

The new two- loop terms and Z^ 2 \ which explicitly depend on the top-mass 
renormalization scheme, can be written as 

r( , = Wr+mn ^ + ^iy r , 

z« = ^ + JVzln g + ( ln J-i) rz . (U, 

The expressions for Wr, N r , Wz and N z are the same in both schemes (up to the 
different m t in the definition of x and y): 

(764 + 3927 y - 9138 y 2 + 6175 y 3 ) y 



W r = 



729 (y-1) 4 
4(32 + 18|/-132 j / 2 + 95i/ 3 )y L , / y-l \ 

8i( y -i) 4 l2 v // ; 

4 (-110 + 797 w- 1233y 2 + 602y 3 + 88?/ 4 )|/ , 

f my 

243 (y - l) 5 

8 (16 + 5y -57y 2 + 54 y 3 )y 2 
^ ln •'/• 



81 (y - 1Y 

4 (263 -486 y + 243 y 2 + 88 y 3 )y 8 (16 + 5y - 57 y 2 + 54y 3 ) y 

-/Vf — ^ — E ■III'?/, 

243 (y-1) 4 81 (y-1) 5 

= 4(-3 + y)*y (2-yW /y-l\ 

3(y-i) 2 (y-i) 2 2 v v ; 
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(2 + 9y-3y 2 )xy A 2xy 2 
•- o my o In 

3(y-l) 3 (I/-1) 3 



where the function Li 2 (z) is defined as 

Li 2 (z) = - f ^ln(l-t). (13) 
jo t 

The expressions for Tr and Tz depend on the renormalization scheme used for m t : 



MS - scheme : T r = 0, 
T z = 0, 



4 (31 - 59 y + 31 y 2 + 9 y 3 ) y 16 (1 - 3 y 2 + 3 y 3 ) y 

27(y-l) 4 9(y-l) J 
(3 -Ay + y 2 )xy 2xy 

(y- Q 3 " (y-i) c 



pole — scheme : 1 r = 7777 774 777 77^ In y, 

t z = ^-;y^ 3 ^y + _^_ lny . ( i 4) 



3.3. Impact of the two-loop contributions on Cio,h 

In this section we briefly illustrate the impact of the two-loop corrections presented 
in this letter on Cio,h(a0- We introduce a rescaled Wilson coefficient (see eq. (4)) 

1 4tt 

In fig. 2 we plot the quantities 

^cSJh(^) and ij (cg H (M + ^^^o!h(m) , (16) 

i.e. two approximations of Cio,h as a function of the charged Higgs boson mass M# 
for the MS- and for the pole mass scheme of the t-quark mass. As input parameters 
we use a s (M z ) = 0.119, m]? olc = 178.0 GeV, M w = 80.4 GeV and s 2 w = 0.231 
[15,16]. The upper frame shows these quantities at the relatively low matching scale 
Hw = Mw- As in this case mf olc and fn t (fJ,w) are numerically almost identical, the one- 
loop approximations (dotted and dashed lines) are close to each other. The inclusion 
of the two-loop corrections, however, considerably lowers the (absolute) size of the 
coefficient for all values of M H considered. In the lower frame a higher matching scale 



of /Jw = 300 GeV is chosen. As in this case m\ oe and m, t (f^w) differ considerably, 
the renormalization scheme dependence of the one-loop results is rather large. When 
taking into account the two-loop corrections (solid and dash-dotted lines), the scheme 
dependence is drastically reduced. 

Looking at the renormalization group equation (RGE) [17] for Cio,h, one finds that 
Cio,h does not run, i.e. 

Cio,h(a*&) = Cio,h(a*w) , (17) 

where the low scale fj, b is of the order of m b . In fig. 3 we show the dependence of 
Cio,h (/-*&) ° n the matching scale nw for Mh = 300 GeV. It can be clearly seen that 
the inclusion of the two- loop contributions drastically lowers the dependence on /i w . 
For /i w > 250 GeV, Cio,R-(/-t&) at two-loop precision is nearly /i^-independent. For /i w 
between M w and 250 GeV the two-loop Wilson coefficient varies about ±4%, whereas 
the corresponding one-loop coefficient varies about ±11%. 

To summarize: In this letter we have presented QCD corrections to the charged Higgs 
induced contributions to the Wilson coefficients Cg(/iw) an d C 10 (iiw) in type-I an d 
type-II 2HDMs. These two-loop results are important ingredients for complete NNLL 
calculations of various observables related to the decay B — > X s l + l~ in these models. 

Just before submitting the present paper, we became aware of the PhD thesis of Ch. 
Bobeth (http://tumbl.biblio.tu-muenchen.de/publ/diss/ph/2003/bobeth.pdf), where 
the two-loop results for the charged Higgs boson contribution to Cg and C w are con- 
tained. We have checked that our results agree. 

We would like to thank K. Bieri and D. Wyler for helpful discussions. S.S would like 
to thank M. Misiak and J. Urban for fruitful discussions and advice regarding the 
technical details of the two-loop calculations. This work is partially supported by: 
the Swiss National Foundation; RTN, BBW-Contract No. 01.0357 and EC-Contract 
HPRN-CT-2002-00311 (EURIDICE). 
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Figure 2: Dependence of the rescaled Wilson coefficient CichO^w) ( see eq. (15)) on 
the charged Higgs boson mass M H at the matching scale nw = M w (upper frame) and 
fj,w — 300 GeV (lower frame). The dashed (dotted) line is the one- loop contribution 
expressed in MS -scheme (pole-mass scheme) of the t-quark mass, while the solid (dash- 
dotted) line includes the two-loop corrections in the respective scheme. 
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Figure 3: Dependence of the rescaled Wilson coefficient Cio,h (/•*&) on the matching scale 
Hw (see eq. (15)) for M H = 300 GeV. The dashed line shows the one-loop contribution 
expressed in MS -scheme for the t-quark mass, while the solid line includes the two-loop 
corrections in the same scheme. 
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